Magnetic field effects on the thermonuclear combustion front of 
Chandrasekhar mass white dwarfs 

Cristian R. Ghezzi 1 , Elisabete M. de Gouveia Dal Pino 1 & Jorge E. Horvath 1 
Received ; accepted 



1 Instituto Astronomico e Geofisico, University of Sao Paulo, Av. Miguel Stefano, 4200, 
Sao Paulo 04301-904, SP, Brasil; E-mail: cghezzi@jet.iagusp.usp.br, 



- 2 - 
ABSTRACT 

The explosion of a type la supernova starts in a white dwarf as a laminar 
deflagration at the center of the star and soon several hydrodynamic instabilities, 
in particular, the Rayleigh- Taylor instability, begin to act. A cellular stationary 
combustion and a turbulent combustion regime are rapidly achieved by the 
flame and maintained up to the end of the so-called flamelet regime when 
the transition to detonation is believed to occur. The burning velocity at 
these regimes is well described by the fractal model of combustion. Using 
a semi-analytic approach, we describe the effect of magnetic fields on the 
fractalization of the front considering a white dwarf with a nearly dipolar 
magnetic field. We find an intrinsic asymmetry on the velocity field that may 
be maintained up to the free expansion phase of the remnant. Considering the 
strongest values inferred for a white dwarf's magnetic fields with strengths up to 
10 8 — 10 9 G at the surface and assuming that the field near the centre is roughly 
10 times greater, asymmetries in the velocity field higher than 10 — 20% are 
produced between the magnetic polar and the equatorial axis of the remnant 
which may be related to the asymmetries found from recent spectropolarimetric 
observations of very young SN la remnants. Dependence of the asymmetry with 
white dwarf composition is also analyzed. 

Subject headings: Thermonuclear combustion: theory, general, fractal 
model, turbulent combustion; - stars: white dwarfs, magnetic fields, 
supenovae, SN la 
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1. Introduction 

The explosion of a type la supernova begins with the combustion, at the center of 
a Chandrasekhar mass white dwarf of carbon-oxygen (C+O) or oxygen-neon-magnesium 
(O+Ne+Mg) fuels. The heat is transported mainly by conduction due to degenerate and 
completely relativistic electrons as a subsonic deflagration wave propagatating outwardly of 
the star. 

The deflagration front born laminar is subject to several hydrodynamic instabilities 
such as Landau-Darrieus (LD) and Rayleigh Taylor (RT) instability (Arnett & Livne 1994, 
Khokhlov 1993) that produce an increment of the area at which the nuclear reactions take 
place. This causes an increase of the nuclear energy generation rate and consequently an 
acceleration of the front. There are two scale ranges that can be distinguished, at the 
smallest or "microscopic" perturbation scales the LD predominates over the RT instability, 
while for the greatest scales 10 5 cm) the RT instability is more important^. The 
combustion front is stabilized by the merging of cells, the formation of cusps, and the 
expansion of the exploding star, which leads to the formation of a cellular structure at 
microscopic scales. LD instabilities lead to an aceleration no higher than 2% near the center 
of the star because they are nonlinearly stabilized (Khokhlov 1995). On the other hand, at 
the macroscopic scales there is a critical wavelength above which the nonlinear stabilization 
failsfl Zel'dovich et al. (1985). 

The bubbles that grow due to RT instability are also subject to Kelvin-Helmoltz (KH) 

2 This depends on the value of the Froude number F = vf am /gL (where L and vi am are 
defined below), then if F < 1 RT instability predominates over LD. 

3 At scales greater than this wavelength the Froude number is F <C 1, and the fluid is 
fully turbulent. 
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or shear instability when nonlinear stabilization fails. The onset of the KH instability marks 
the transition to the fully developed turbulence regime. During this, fluid motions are 
characterized by the formation of a turbulent cascade in the inertial scales where viscous 
dissipation is not important. This turbulence can be described by the Kolmogorov's scaling 
law. 

The acceleration of the turbulent thermonuclear flame due to the action of the several 
hydrodynamics instabilities can be described by the fractal model first introduced by 
Woosley (1990), Timmes & Woosley (1992), Niemeyer & Hillebrandt (1995) and Niemeyer 
& Woosley (1997). The idea of applying fractals to the combustion fronts is useful because 
the statistic properties of the surface change self-similarly on the spatial scales. The 
wrinkled surface S behaves like a fractal S oc R D , where D is the fractal dimension of 
the surface, with 2 < D < 3, and R is the mean radius of the wrinkled surface (Filyand, 
Sivashinsky Sz Frankel 1994). There is a scale range at which the similarity holds, called 
the similarity range, that goes from l m in to L. The increase of the area at which the nuclear 
reactions occur induces an increase in the propagation velocity that also behaves like a 
fractal^ (see Niemeyer & Hillebrandt, 1995). This fractal model is supported by laboratory 
experiments involving different gas mixtures (Gostintsev, Istratov & Shulenin 1988). 

The propagation of the front subject to LD instability is well described by the 
Sivashinsky equation (Sivashinsky 1977, 1983) that also includes thermodiffusive and 
acceleration effects, or by its gereralization (Frankel 1990). Numerical studies of wrinkled 

4 We refer here to the fractal dimension (and not to the Hausdorff-Besicovitch dimension) 
of an attractor as it is commonly defined in astrophysical works (Barnsley 1988). The 
wrinkled surface is not a fractal in the mathematical sense, it is converging to its attractor 
asymptotically with time, but it is possible to treat the surface as a fractal as long as 
I mini L <C 1, for ex., l min /L ~ 10~ 5 at p ~ 10 8 gem -3 (see below). 
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surfaces utilizing these equations are able to reproduce the front on a wider range of spatial 
scales than direct hydrodynamic simulation and allow the study of the fractal properties of 
the flame (Blinnikov & Sasorov 1996). These simulations of the Sivashinsky equation show 
that the fractal growth of the front increases its velocity (Filyand, Sivashinsky & Frankel 
1994). 

Several physical effects are yet to be considered for a complete solution of the SN la 
problem. In particular, the effects of magnetic fields which are known to be present in SN 
la progenitors (believed to be white dwarfs in binary systems) have not been considered. 
We will incorporate here the effects of magnetic fields on the fractal growth of the turbulent 
combustion front and show that they can break the spherical symmetry of the explosion. 

2. Formulation of the problem 

Ginsburg (1964) and Woltjer (1964) proposed that the magnetic flux of a star is 
conserved during its evolution and collapse, so that strong magnetic fields would be 
generated in degenerate stars. Hence, a main sequence star with R ~ 10 10 cm and 
B ~ 10 - 10 5 G would collapse to form a white dwarf with R ~ 10 8 cm and B ~ 10 5 - 10 9 G 
on its surface. Inferred magnetic field strengths of known magnetic white dwarfs range from 
~ 10 6 to 10 9 G (Jordan 1992). There is also an important percentage of known white dwarfs 
in cataclysmic variables with magnetic fields ranging from ~ 10 6 to 10 8 G (see Chanmugan 
1992, Liebert and Stockman 1983). There are several models of white dwarfs with strong 
magnetic fields (Jordan 1992, Putney & Jordan 1995, Martin & Wickramasinghe 1984), 
which assume centered dipolar, off-centered dipolar, quadrupolar or dipolar+quadrupole 
magnetic field configurations. Although the origin of the fields in these stars remains 
unclear it is frequently assumed that the fields are primordial. In fact, Wendell, Van Horn 
& Sargent (1987) carried out detailed calculation of the time evolution of a white dwarf's 
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magnetic field and found the decay time of the field is always longer than the typical cooling 
times of white dwarfs (r ~ 10 10 yr). We will assume in this letter that the progenitor star 
of an SN la is a magnetic white dwarf with a centered dipolar magnetic field. Considering 
magnetic fields with the strongest inferred values, i.e., in the range 10 8 — 10 9 G, and using 
the model of Wendell, Van Horn & Sargent (1987), we find magnetic field strengths near 
the center of the star which are in the range ~ 10 9 — 10 10 G. We thus must expect that 
when the explosion begins the flame will propagate parallel to the magnetic fields lines at 
the magnetic poles of the star, and perpendicularly to the field lines at the equator. 

An estimate of the effective turbulent speed u t of the flame was obtained by Damkohler 
(1939), who proposed that u t ~ v{L), where v(L) is the average of the turbulent velocity 
fluctuations on the largest hydrodynamical scale L. Later Karlovitz, Denniston, & Wells 
(1951) derived a statistical approach based on turbulent velocity correlations that fits better 
with laboratory experiments. The formula simplifies to u t = \2ui arn y(L)] l l 2 in the limit 
v(L) vi am , where vi am is the laminar velocity. Niemeyer & Hillebrandt (1995) generalized 
those previous models and found 



where the exponent n is arbitrary, with n = 1 for Damkohler's model and n — 1/2 for 
the model of Karlovitz et al.. Although numerical simulations of thermonuclear flames 
are in good agreement with the Karlovitz et al. model (see Khokhlov 1993), Niemeyer & 
Hillebrandt use their own version to derive the scaling law for the turbulent flame speed. 
We will give here a slightly different derivation of the scaling law for u t . The turbulent 
motions of the fluid can be described by the Kolmogorov's scaling law for incompressible 
turbulent velocity fluctuations which gives (Landau & Lifshitz 1959) 



this scaling is valid in the inertial range (77 <C I <C L) between the dissipation scale r] 



Ut(l) = Vl am [v{l)/Vi am \ 



n 



(1) 




(2) 
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and the length scale L at which the turbulent velocities are in equilibrium with the RT 
instability. So if we assume that, at the largest hydro dynamic scales which are subject to 
the RT instability, the turbulent velocity is in equilibrium with the velocity of RT bubbles 
vrt(L) ~ v{L), where vrt(L) oc (gL) 1//2 , we can use eq. (2) to obtain the velocity of the 
turbulent fluctuations at the lower scales v(l) = f^ r (/)(L//) 1 / 6 . From this it is clear that 
v(l) > vrt(1) for / < L. Therefore, the turbulent cascade will dominate on all scales below 
the largest RT unstable wavelength (see also Niemeyer & Hillebrandt 1995). Using the last 
formula and eq. (1), we obtain u t (L) = Ui am [vRT{l)/u lam ] n {L/l) n / & . Following Timmes k 
Woosley (1992), we will use as a lower cutoff the minimum length scale l m i n that can deform 
the laminar flame front. This length is given by the condition v RT {l min ) = V[ am . Therefore 
the scaling law for the turbulent flame speed is 



In a fractal description of the turbulent flame propagation, eq. (3) gives the fractal 
velocity of the front Vf rac = vi am (L j l m in) D ~ 2 (Niemeyer & Woosley 1997, Niemeyer & 
Hillbrandt 1995). Comparing this with eq. (3), we see that the fractal excess is D — 2 = 0.5 
for n = 1, and 0.25 for n = 1/2, which is in agreement with the fractal dimensions inferred 
from numerical studies of SN la explosions (Blinnikov, Sasorov & Woosley 1995). 

In the equation above we must determine l min and L. As we have seen, the characteristic 
RT velocity of the growing modes must be larger than or equal to the laminar deflagration 
velocity vi am , or v RT {l m in) = vi am . This implies that l m inn RT {l m in) = vi am , where n RT {l) is 
the inverse of the characteristic RT time n R T = (l/2n)J gkAp/2p. This results 



where g is the acceleration of the gravity, Ap = p u — Pb is the difference between the 
densities p u and pb of the unburned and burned fuel, respectively, and p = (p u + Pb)f2. We 




n/2 



(3) 



mm 




(4) 
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denote / as l vo \ hereafter, where the subscript "poZ" stands for polar since it gives the value 
of l min at the magnetic poles of the star where the propagation of the flame is parallel to 
the magnetic field lines, and therefore unaffected by them. 

The value of L has been determined from the minimum between the maximum spatial 
extent of the density inversion and the radius of the flame Rf. The density inversion is due 
to electron capture in an isobaric nuclear statistical equilibrium environment that causes 
the number of electrons to decrease and the density to increase. We here use L = Rf 
because, for densities <; 10 9 gem -3 , the density inversion scale is greater than the stellar 
radius (Timmes & Woosley 1992). The highest unstable mode L has the highest growing 
time since £_rt(0 = (47r//(?) 1 / 2 (p/Ap) 1 / 2 < t RT (L), for l po i < I < L. For example, near 
the center of the star L ~ 3 x 10 7 cm and the perturbation takes t RT {L) ~ 0.5 s to grow 
to an amplitude of the order of its wavelength. On the other hand, for l pol ~ 10 5 cm 
tRrilpoi) = 0.03 s and the perturbation grows faster. Thus with a deflagration velocity 
viam ~ 10 7 cms _1 , the flame must travel only ~5x 10 6 cm for the slowest perturbations, 
and ~ 3 x 10 5 cm for the fastest ones in order to the perturbations to grow to amplitudes 
of the order of their lengths. This suggest that all the perturbation scales grow very fast in 
the combustion front and therefore, they are all important. 

Our analysis is applicable only for densities p > 10 7 gcm~ 3 , since the turbulent motions 
will destroy the corrugated flamelet regime for lower densities (Niemeyer & Woosley 1997, 
Hillebrandt & Reinecke 2000). 

3. Asymmetric velocity field 

Near the center of the star the magnetic pressure (B 2 /8tt ~ 10 18 ergcm~ 3 ) is much 
smaller than the gas pressure behind the front (p ~ 10 27 erg cm^ 3 ). Therefore, it is 
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reasonable to assume that the scaling law for the flame speed (eq. 3) will not be modified by 
the presence of the magnetic field. This means that the fractal dimension can be considered 
independent of the magnetic field. On the other hand, the presence of magnetic fields will 
shift the lower unstable scales at the equator of the star with respect to the poles, as we will 
see, thus producing different fractal velocities at the magnetic poles and the equator of the 
star. We note also that the turbulent fluid motions can only scramble the field inside the 
"flame brush", that is, inside the RT foam that mixes unburned and burned material, or 
behind the front, and therefore, they will not affect the large scale magnetic field geometry 
of the progenitor. 

The dispersion relation for the RT instability of the front with a magnetic field 
transverse to the direction of the gravity and to the propagation of the flame, in a 
heterogeneous inviscid plasma of zero resistivity is (Chandrasekhar 1961): 



nRTB = iv k {^-^) (5) 

here k = 2ir/l, and B is the magnetic field strength. It is, therefore, possible to define an 
effective surface tension, T e ff = B/(2nk), for which the magnetic field will have a stabilizing 
effect over perturbations with wavelengths below ~ B 2 /(gAp). However, if the velocity of 
the flame is taken into account another cutt-off length appears, l eq , which is given by the 
condition l eq nRTB(hq) = v iam- From this equation the minimum scale of the self-similar 
range in presence of magnetic field is 

, ._ Mg + Warn) ((]] 

Leq ~ gAp W 

here eq stands for "equator" because l eq gives the minimum instability scale at the magnetic 
equator of the star where the propagation of the front is perpendicular to the field lines. 
Inserting eq. (6) in the definition of the fractal velocity (eq. 3) we have 

leq - 



Veq = Vl am ( — ) (7) 
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for the fractal velocity in the magnetic equator of the star, while, as we have seen before, 
eq. (3) gives v po i = u t for the magnetic poles, 

/ l \ D-2 

Vpol = Vl am — (8) 
V l-pol / 

with l po i given by eq. (4). 

Taking the ratio between the velocities v po i and v eq , we obtain 

Vpol ( B 2 /8ir + pv? am /2 \ D -i 

v eq \ pvf am /2 J ■ W 
The relevant data for do the calculations can be found in Timmes & Woosley (1992). 
We will use a fractal dimension D = 2.5 as in Blinnikov, Sasorov & Woosley (1995) and 
references therein. 

Figure 1 displays the percentage of asymmetry in the velocity field for two white 
dwarfs with different compositions: (a) X( l2 C) = 0.5 X( 16 0) = 0.5 , Ap/p = 0.426, 
v lam = 2.33 x lO^ms" 1 and (b) X( 12 C) = 0.2 X( 16 0) = 0.8, Ap/p = 0.415, 
viam — 0.415 x 10 5 cms _1 , and fuel densities p ~ 10 8 gcm -3 which are encountered by the 
front at a middle radius distance ~8 x 10 7 cm from the center of the star. We see that the 
asymmetry is very sensitive to the composition of the progenitor and is quite insignificant 
for type (a) progenitor. If at these densities, magnetic fields as high as ^ 6 x 10 8 G are 
present at this radius, the velocity field has an asymmetry ^ 10 % for progenitors of type 
(b). With surface magnetic fields of the order of 10 9 G, and fields ^ 6 x 10 9 G at the 
interior, these progenitors could suffer huge asymmetries ^ 120% (see Fig. 1). Progenitors 
with compositions 16 20 Ne show even higher asymmetries (see Ghezzi, de Gouveia Dal 
Pino & Horvath 2001) because the laminar velocity in 16 O 20 Ne is smaller than in 12 C 16 0. 
Also, we note that for regions very close to the center of the star (at higher densities) the 
asymmetry percentages are smaller than 5 %. 

When the flame reaches the lates stages of the flamelet regime, at a density of 
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~ 5 x 10 7 gcm~ 3 , the magnetic field intensity drops to ~ 10 8 — 10 9 G, however our 
calculations indicate that the percentage of asymmetry is maintained or increased (Ghezzi, 
de Gouveia Dal Pino & Horvath 2000), leading to an asymmetry higher than > 10 % for a 
progenitor with X( 12 C) = 0.2 X( 16 0) = 0.8 with fields higher than 10 s G at the surface of 
the star. 

4. Conclusions and Discussion 

We have found that an asymmetric field velocity caused by the presence of a dipolar 
magnetic field during the fractal growth of the deflagration front of a type la supernova 
can lead to the formation of a prolate remnant. The magnetic field introduces an effective 
surface tension in the equator of the white dwarf progenitor that reduces the velocity of the 
combustion front, v eq , with respect to the velocity at the poles, v po i, so that v vo \ > v eq by a 
percentage of 10 % to 20 % for progenitors with a composition X( 12 C) = 0.2 X( 16 0) = 0.8 , 
Ap/p = 0.415, and surface fields ~ 10 8 G (type (b), Fig. 1). This leads to prolate explosions 
along the magnetic poles. Lower magnetic field strengths have no detectable effects on the 
explosion. As only a small fraction of the observed white dwarfs must have magnetic fields 
as high as 10 8 G, asymmetric explosions are not expected to occur very frequently. If there 
is no transition to detonation the asymmetry will be probably maintained during the free 
expansion phase of the supernova remnant, because the expansion velocity is constant at 
this phase. Detonation in an SN la is still controversial, but if a transition to detonation 
occurs the remnant can be symmetrized very fast (Ghezzi, de Gouveia Dal Pino & Horvath 
2001). 

Recent spectropolarimetric observations have revealed a linear polarization component 
in the radiation of very young SN la remnants, which suggests that prolate atmospheres 
with asymmetries ~ 20 % are producing it (see Leonard, Filippenko, & Matheson 1999, 
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Wang , Wheeler & H officii 1997). The model presented here offers a plausible explanation 
for such observations and provides new motivation for theoretical studies of supernovae 
involving magnetic fields. However, its predictions must still be confirmed through fully 
numerical simulations of the explosion in the presence of magnetic fields. 

This paper has benefited from many valuable comments by the referee F. X. Timmes, 
and discussions with M. Diaz. C.R.G., E.M.G.D.P. and J.E.H. have been partially 
supported by grants of the Brazilian Agencies FAPESP and CNPq. 
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Fig. 1. — Asymmetry percentage in the field velocity for two progenitors with initial 
compositions X( 12 C) = 0.5 X( 16 0) = 0.5 and X( 12 C) = 0.2 X( 16 0) = 0.8 , at p ~ 10 8 
g cm~ 3 , as a function of the magnetic field strength at a radius ~ 8 x 10 7 cm. The asymmetry 
percentage is given by: 100 x (^v po i/v eq — l) where v po i/v eq is obtained by eq. (9). 



